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Abstract 

We study density correlation functions for an impenetrable Bose gas in a finite box, with Neu- 
mann or Dirichlet boundary conditions in the ground state. We derive the Fredholm minor determi- 
nant formulas for the correlation functions. In the thermodynamic limit, we express the correlation 
functions in terms of solutions of non-linear differential equations which were introduced by Jimbo, 
Miwa, Mori and Sato as a generalization of the fifth Painleve equations. 



Research Fellow of the Japan Society for the Promotion of Science. 



1 Introduction 

In the standard treatment of quantum integrable systems, one starts with a finite box and imposes 
periodic boundary conditions, in order to ensure integrability. Recently, there has been increasing 
interest in exploring other possible boundary conditions compatible with integrability. 

With non-periodic boundary conditions, the works on the Ising model are among the earliest. By 
combinatorial arguments, McCoy and Wu [jlj studied the two-dimensional Ising model with a general 
boundary. They calculated the spin-spin correlation functions of two spins in the boundary row. 
Using fermions, Bariev § studied the two-dimensional Ising model with a Dirichlet boundary. He 
calculated the local magnetization and derived the third Painleve differential equations in the scaling 
limit. Bariev || generalized his calculation to a general boundary case. In the Neumann boundary 
case, he also derived the third Painleve differential equations in the scaling limit. Sklyanin Q began 
a systematic approach to open boundary problems, so-called open boundary Bethe Ansatz. Jimbo 
et al H calculated correlation functions of local operators for antiferromagnetic XXZ chains with a 
general boundary, using Sklyanin's algebraic framework and the representation theory of quantum 
affine algebras. 

Sklyanin 0] explained the integrability of the open boundary impenetrable bose gas model, using 
boundary Yang Baxter equations. In this paper, we will study field correlation functions (density 
matrix) for an impenetrable bose gas with Neumann or Dirichlet boundary conditions. Schultz || 
studied field correlation functions for an impenetrable bose gas with periodic boundary conditions. 
He discretized the second quantized Hamiltonian and found that the discretized Hamiltonian was 
the isotropic XY model Hamiltonian. He diagonalized the discretized Hamiltonian by introducing 
fermion operators. Using the N particle ground state eigenvector for the discretized Hamiltonian, 
Schultz derived an explicit formula of correlation functions for an impenetrable bose gas in the 
continuum limit. Lenard Q] pointed out that Schultz's formula could be written by Fredholm minor 
determinants. Therefore this formula is called Schultz-Lenard formula. In this paper, we will derive 
Schultz-Lenard type formula for Neumann or Dirichlet boundary condition. Following Schltz, we 
employ two devices. We consider the N particle ground state of the discretized Hamiltonian. We 
then fermionize the discretized N particle system by using the Jordan- Wigner transformation. In 
the continuum limit, we derive the Fredholm minor determinant formula for correlation function, 
which has the integral kernel : 
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[L : box size, N : the number of particles, e = + : Neumann, e = — : Dirichlet.) 

Jimbo, Miwa, Mori and Sato developed the deformation theory for Fredholm integral equation 
of the second kind with the special kernel sm j^~f ^ . They introduced a system of nonlinear partial 
differential equation, which becomes the fifth Painleve in the simplest case. They showed that 
the correlation functions without boundaries was the r-function of their generalization of the fifth 
Painleve equations. In this paper, we express the correlation functions for Neumann or Dirichlet 
boundaries in terms of solutions of Jimbo, Miwa, Mori and Sato's generalization of the fifth Painleve 
equations, hereafter refered to as the JMMS equations. In the thermodynamic limit (N, L — > oo, ^ : 
fixed), we reduce the differential equations for correlation functions with Neumann or Dirichlet 
boundaries to that without boundaries, using the reflection relation between two integral kernels 
sm(x-x ) sm(x+x ) ^ sm(x-x ) rj,^ point correlation function with Neumann boundary is 

x—x' x+x' x—x' r J 



described by the equation ( |2.29| ) and ( 2.3C| ). In the case with boundary, the differential equation 



for the two point correlation function cannot be described by an ordinary differential equation. We 
need three variable case of the JMMS equations. 

Physically, the long distance asymptotics of the correlation function are interesting. The long 
distance asymptotics of the ordinary differential Painleve V is known. But, for many variable case, 
the asymptotics of the JMMS equations are not known. Therefore we cannot describe the long 
distance asymptotics of the correlation functions with boundary in this paper. To evaluate the 
asymptotics of the solution of the JMMS equation is our future problem. 

Now a few words about the organization of the paper. In section 2, we state the problem and 
summarize the main results. In section 3, we derive an explicit formula for the correlation functions 
in a finite box. In section 4, we write down the differential equations for the correlation functions 
in the thermodynamic limit. 



2 Formulation and results 

The purpose of this section is to formulate the problem and summarize the main results. The 
quantum mechanics problem we shall study is defined by the following four conditions. Let N € 
N,(N > 2),L € R,$o,$ L G R- 

1. The wave function ipN,L = ?Pn,l(xi, • • • , xnWo, $l) satisfies the free-particle Schrodinger equa- 
tion for the motion of N particles in one dimension (0 < (xi ^ Xj) < L). Here the variables 
x\, • ■ ■ , xn stand for the coordinates of the particles. 
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2. The wave function ipN,L is symmetric with respect to the coordinates. 

iPn,l(xi, - ■ ■ ,xn\"&o,&l) = it>N,L(x a (i), ■ • • ,z ct (jv)|#o,$l), ( a G s n)- (2-1) 

3. The wave function satisfies the open boundary conditions in a box < Xj < L, (j = !,■■■ , N) 



7^ #0 j 1pN,L{xi, ■ ■ ■ ,Xn\&0,&l) 



d 

d Xj 



+ 1?L 1pN,L(Xl, ■ ■ ■ ,X N \$Q,$ L ) 



= 0, (j = l,---,iV), 
= 0, (j = !,-••, TV). 



(2.2) 
(2.3) 



4. The wave function tpN,L vanishes whenever two particle coordinates coincide. 

iPn,l(xi, ■ ■ ■ ,xu ■ ■ ■ ,Xj, ■ ■ ■ ,xjv|i?o,i?l) = 0, for Xi = xj. 

In this paper we shall be concerned with the ground state. The wave function is given by 
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(2.4) 



iPn,l(xi, ■ ■ • ,xat|^ ,^l) 



sJv N , L ($ ,d L ) 



det (Xj cos (XjXk) + #o sin (AjXfc)) 



(2.5) 



Here the momenta < Ai < • • • < \n are determined from the boundary condition for ipN,L which 
amounts to the equations 

2LA, + e^Xj) + 6, dL (A,-) = 2irj, (j = 1, • • • , N), (2.6) 

where we set d (\) = Hog (f£M^). We take the branch -vr < 9 d (X) < vr, (d > 0). Here Vn^o^l) 
is a normalization factor defined by 



det f V (A j + £ ^o)(A fe + e'tfo) / L e^-^^dy I 



(2.7) 



The wave function is not translationally invariant and has the normalization, 



/ ••• / dy 1 ---dy N ip N>L (yi,---,y N \'&o,tf L ) 2 = 1. 
Jo Jo 



The following equation holds for any parameter A, 

d 



dx 



■&o (Acos(Ax) + -&o sin(Ax)) 



= 0. 



x=0 



The following equivalent relation holds, 



— + i? L ) (A cos(Ax) + -do sin(Ax)) 



& e 



2iL\ 



x=L 



(A + ^ l )(A + ^q) 
(A-m? l )(A-«V> 



(2.8) 



(2.9) 



(2.10) 



/,From fl2.9|) , ( 2.1 0| ) and Girardeau's observation on fermions and impenetrable bosons correspon- 



dence in one dimension ||, we can show that the wave function ip^L satisfies the above four condi- 
tions. We shall be interested in the field correlation functions (density matrix) given by 

pn,N,L{x\, ■ ■ ■ ,x n \x[, ■ ■ ■ ,x' n \i} ,$ L ) 
(n + N)\ f L rL 



o Jo 



dy n +i ■ ■ ■ dy n+ Ntp n +N,L(xi, ■ ■ ■ ,x n ,y n+ i, ■ ■ ■ ,y n+N \'&o,^L) 



X i>n+N,L(x'i, ■ ■ ■ ,x'n>yn+l, ■ ■ ■ ,yn+N\$0,fiL)- (2.11) 



In this paper, following ||], we reduce our problem to that of discrete M intervals. Set e = M L +1 . 
Let \vi) = I ] , \v2) = ( ) be the standard basis of V = C 2 . Let (vi\, (i = 1,2) be the dual 

W v 1 / 

basis given by (vi\vj) = 5{j, = 1,2). The action of O G End(C 2 ) on (vi\, (i = 1,2) is defined 
by (( Vi \0) \vj) := ( Vi \ (0\v d )) , (j = 1,2). Set \Sl ) = \v^ M and (J2 | = ( Vl \® M . Set 



/ 1 \ ,10 

10/ \00/ \ -1 



(2.12) 



Following the usual convention, we let ^,^,<r| signify the operators acting on the j— th tensor 
component of V® M . Introduce fermion operators -0m j VVt by the Jordan- Wigner transformation 

i>m = 0*1 •••°'m-l ( Am ) i>m = of ■ ■ ■ &m-l<l>m, i™, = 1, • • • , M). (2.13) 

The fermion operators have the anti-commutation relations 

{^+,V>n} = 5 m ,n, {^m, = } = 0- (2.14) 

Here we use the notation {a, 6} = ab + 6a. Set 



n~ M 

\jjf\ X, V , JV,i(emi,---,emjv|'?o J ^L)^mi ■•• < / > m Jv l^o) (2.15) 



mi,---,r?ijv=l 

AT A/ 



IN\ V n ,l(#o,0l) "lm" 

n - M 

(O^mC^O,^)! = \/T7T X 0AT,L(emi,---,emAr|^o,^L)(^o|'/ > mi ■■■<Pm N (2.16) 



II X (A j cos(em i A i ) - tf sin(emjAj)) ^ • • • V4j fi o), 



iV! 

mi,---,mjv=l 



, N M 

-, TT XI (Aj cos(em J A J ) - i? sin(emjAj)) (fiolVWi ••■V , m JV - 
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Using the above vectors, we can calculate correlation functions in the continuum limit as 

Pu,N,l{xi, ■ ■ ■ ,x n \x[, ■ ■ ■ ,X , n \ti ,'&L) 
( L\ N 

= ^^{jfj (^n+NM^O^LMs^Si • • • <i>s n 4>Z4>ti ■■■ ( t ) tj n n+N,M(' & 0,^L)), (2.17) 

where we take the limit M — » oo in such a way that esj — > Xj , etj — > x ^ , (L : fixed) . The equation 
(pT7| ) follows from ( gig ) and (gl9| ). 



iV M 

m j^ t i''"' in ' s i''"' Sn 

x ^n+JV,i(e*i • • • e*n, ewi • • • emjvl^o, &l) ip n +N,L(eh ■ ■ ■ et n , em\ ■ ■ • emjv|i?o, "&£) 

N+n M N+n M 

= II 5Z II H^ol^i •••0m n+Jv ^ 1 ■■■^a n 0ti ■ ■ ■ 4>t n (f>h ■ ■ ■ (t>l n+N \tt ) 
i=l m,i=\ j=l lj=l 

x ipn+N,L{emi ■ ■ ■ em n+N \'d ,'d L ) ip n +N,L{eh • • • e/ n+A r|t? , ^l), (2.18) 

(J2o|0mi • ■ ■ <?W^mi ' ' ' ^mjv l^o) = 1- (2.19) 



This formula (2.17) is our standing point. The case ^Oj^l = corresponds to Neumann boundary 
condition and the case ^Oj^z, = °° to Dirichlet boundary condition. In the sequel, we use the 
following abbreviations. 

Pn,N,L{x\, ■ ■ ■ ,X n \x[, ■ ■ ■ ,x' n \ + ) = pn,N,L( X li'--,X n \x' 1 ,--- ,x' n \0,Q), (2.20) 

Pu,n,l(xi, ■ ■ ■ ,x n \x[, ■ ■ ■ ,x' n \-) = Pti,n,l(xi,- ■ ■ ,x n \x' 1: - ■ • ,o4|oo,oo). (2-21) 

In the sequel, for simplicity, we consider two important case : Neumann boundary conditions and 
Dirichlet boundary conditions. 

Remark. There exists the simple relations between Neumann or Dirichlet boundaries and 
periodic boundaries. We can embed the differential equations for n point correlation functions of 
Neumann or Dirichlet boundaries, to the one for In or 2n — 1 point correlation functions without 
boundaries. 

In section 3, we derive the following formula. 

Theorem 2.1 The correlation functions for an impenetrable hose gas with Neumann or Dirichlet 
boundaries are given by the following formulas. 

( l\ n 

p n ,N,L( x ir--,x' n \x'{,---,x^\e) = (--) J[ sgn(4 - x'j) sgn(4' - x") 

l<j<k<n 
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x det 1 K, 



IT 



e,N,I p 



x'i, x 2 , 
x l5 x 2 , 



x' 
x" 



(2.22) 



where e = ± and < x'^x" < L, (j = 1, ■ ■ ■ , n). i/ere I p zs i/te union of n intervals I p = [x±, x 2 ] U 
■ ■ ■ U [x2n-i) x 2n], where < xi < • • • < x 2ra < L is the re-ordering of x^, ■ ■ ■ , x' n , x'{, • • • , x'^. The 



symbol det 1 - XK £jN j p 



ry^ fY*^ ■ • ■ Of* ^ 



x l5 x 2 , 



j x n 



denotes the n-th Fredholm minor corresponding 
to the following Fredholm type integral equation of the second kind. 

((l-XK £tN , Ip )f)(x)=g(x), (xel p ). 

Here the integral operator K £ ,N,I P is defined by 



(2.23) 



l p 2L [ sin 27;7r(x - y) 



sin ^tt(x + y) 



Using the above Fredholm minor formulas, we can take the thermodynamic limit for correlation 



functions, i.e., N, L — > oo, = p$ : fixed. 

Corollary 2.2 The correlation functions for an impenetrable bose gas with Neumann or Dirichlet 
boundaries are given by the following formulas in the thermodynamic limit. 



Pn( x 1) ' ' ■ j ^nl-^l ) ' ' ' > ' c nl £ ) 



lim 

N,L— >oo, 77=po 



Pn,N,L(x±, " " " , x 'n\ x l-> ' ' ' ' x n\ £ ) 



(2.25) 



l<J<fc<" 



T* rv» ■ ■ ■ 

t" -r" ... -r" 

1 ? J/ 2) i^n 



where < x'^x'- < +oo, (j = l,---,n). The symbol det 1 — XK £ j p 
represents the n-th Fredholm minor corresponding to the following Fredholm type integral equation 



i x r. 



ry.ll 



of the second kind, 

((l-\k eJp )f)(x)=g(x), (xel p ), 
where the integral operator K £ ,l p is defined by 



(Ksjj) (*) = / 



smp Q ir(x - y) sinp 7r(x + y) 



x-y 



+ e- 



x + 



f(y) dy. 



(2.26) 



(2.27) 
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In the sequel, we choose such a scale that irpo = 1. 

In section 4, we derive the differential equations for the correlation functions. Jimbo, Miwa, 
Mori and Sato || introduced the generalization of the fifth Painleve equations, hereafter refered to 
as the JMMS equations. Their simplest case is exactly the fifth Painleve equation. We reduce the 
differential equations for Neumann or Dirichlet boundary case to that for without-boundary case, 



using the reflection relation in lemma 4.2. For n = 1 and Dirichlet boundary case : 

^(0^1-) = 0. (2.28) 

Next we explain n = 1 and Neumann boundary case. The differential equation for /3i(0|x|+) is 
described by the solutions of the Hamiltonian equations which was introduced in || as the special 
case of the generalization of the fifth Painleve equations. We cannot describe the correlation function 
pi(0|x|+) in terms of the fifth Painleve ordinary differential equation. We need the many variable 
case of the JMMS equations. 

— lo gjOl (0|x|+) = H 2 (-x,0,x). (2.29) 
ax 

Here -ff 2 (ao, a ii a 2) is the coefficient of the following Hamiltonian 

H = H (ao 1 ai 1 a2)da + Hi(a ,ai,a2)dai + H 2 (a ,ai,a2)da2 
= ~ ^( r +i r -i ~ r +l r -j)(r+jr-i - r + if^ j )d\og{a j - a x ) 

j=0,2 

- (r +0 f^2 ~ f+2r-o)(f +0 r-2 ~ r + 2f- )dlog(a - a 2 ) (2.30) 



+ ir + ir-idai + i ^ { r +j^-j ~ ^j r -j)daj - dlog(a - a 2 ). 

j=0,2 

Here the functions r ±j = r ±J (a , «i, a 2 ), (j = 0,1,2), f ±0 = f ±0 (a , at, a 2 ), f ±2 = r ±2 (a , a x , a 2 ) 
satisfy the Hamiltonian equations 

dr ±j = {r ±j ,H}, (j = 0, 1, 2), df ±0 = {f ±0) H}, dr± 2 = {f ±2 ,H}. (2.31) 

where the Poisson bracket is defined by 

{r+i,r_i} = 1, {r+o^-o} = {f +0 ,r_ } = 1, {r +2 ,f_ 2 } = {f +2 ,r_ 2 } = 1. (2.32) 

This Hamiltonian H depends on odd number variables ao, ai, a 2 . In the case without boundary, the 
differential equations for the correlation functions are described by the Hamiltonian equations which 
depend on even number of variables |8|. Therefore this point is new for Neumann boundary case. 

In the general case, we can embed the differential equations for n point correlation functions of 
Neumann or Dirichlet boundaries, to the one for 2n or 2n — 1 point correlation functions without 
boundaries. 
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Theorem 2.3 In the thermodynamic limit, the differential equation for the correlation functions 
becomes the following. 

dlogp n (x[,-..,x'Jxl---,x^\e) = (l-n)u £jp ^+J2 £ 4^ (£) > ( 2 - 33 ) 

where we denote by d the exterior differentiation with respect to x\, ■ ■ ■ , x' n , x'{, ■ ■ ■ , x'^. Here the dif- 
ferential forms uj £ j p (X) and uff ^(A) are defined in Proposition and Proposition 4-4, respectively. 



XTie differential forms cu £ j p (X) and lo k £ f k (A) are described in terms of solutions of the generalized 
fifth Painleve equations which were introduced by M. Jimbo, T. Miwa, Y. Mori and M. Sato J^/. 
Both Neumann and Dirichlet boundary conditions, lo £ j p (X) and OJ^f \X) are described by the same 
solutions of the same differential equations. 



Physically, it is interesting to derive the long distance asymptotics of correlation functions : 

p n (x' 1 ,---,x' n \x'(,---,x^\e). (2.34) 

^From the above theorem, we can reduce the evaluation of the asymptotics to the following two step 
problem. 

1 . Evaluate the asymptotics of the solution of the generalized fifth Painleve introduced in ||| . (For 
our purpose, we only have to consider the special solution related to the correlation functions 
for the impenetrable Bose gas without boundary. ) 



2. Determine the asymptotic solutions of the differential equation ( 2.33 ) under the appropriate 



initial condition. (The main point is to determine the constant multiple in the asymptotics. ) 

In the case reducible to an ordinary differential equation, the above two problems have been 
already solved. Jimbo, Miwa, Mori and Sato || considered the problem 1 of the correlation functions 
for the impenetrable Bose gas without boundary. McCoy and Tang [|l0| generalized the asymptotic 
formulas M to the 2-parameter solution of Painleve V, which is analytic at the origin. Vaidya and 



Tracy [11|, [12| considered the problem 2 of two-point correlation functions for the impenetrable 
Bose gas without boundary. (The pioneering work for Ising model was done by McCoy, Tracy and 
Wu fi~3|] .) In our case, to evaluate the asymptotics of pi(0|x|+), we have to consider the case of three- 
variables. However the asymptotics in many variable case is a non-trivial open problem. Therefore 
the above two problems for many variables case are our future problems. 
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3 Predholm minor determinant formulas 



The purpose of this section is to give a proof of Theorem 2.1. Set V = C 2 . For s = ±, define 
operators ri + (6,e),r](6,e) acting on V® M by 

M 

= E (e- ime + ee ime )v4, (3-1) 

m=l 

M 

n(6,s) = (e im0 + ee- imd )i; m . (3.2) 

m=l 

In the sequel we use the notation O^m = m+i 71 "- ^he operators t/ + (6> Mi a/, e), rj^O^MjS) have the 
following anti-commutation relations for e = ±, — M < /i < M. 

{v + (,0iJ<,M> £ )i r l + (9v,M,£)} = {v (0n,M,£),ri (Ov,M,e)} = 0, (3.3) 
|?? + (^,M,e),?? (0 v ,Af,e)} = 2(M+l)((5 M)I/ + e* M) _„). (3.4) 

In the sequel we use the following abbreviations. 

|Ojv,m(+)) = |njv,Af(0,0)) 5 |Ojv,m(-)} = |njv,Af(oo,oo)), (3.5) 
<0iV,M(+)| = <fijv,M(0,0)|, (n NtM (-)\ = {n NM (oo,oo)\. (3.6) 

Using the operators r] + (8,e), rj(6,e), we can write 



1 



n N)M (s)) = Jj^^v + (Oi,M,e)v + (02,M,e)---7 1 + (e NM ,emo), (3.7) 



(SIn,m( £ )\ = (2eL) N ^oMfljv.Af , e) ' ' ' v(Q%M,z)n{6i,Mi £ )- ( 3 -8) 

The operators i) + (9 ft M,£), f]{Qfx,Mi £ ) ac t on t ne vectors in^M^)), (^jv,m( £ )| as follows. 

For < JV, 77 + (^,M,e)|O iVi M(e)} = 0, (fijv,Af (e)|»7 (^Af,e) = 0. (3.9) 
For > N, (n NM (e)\v + {0^,M^) = ^ V {0„,M> £ )\ n N,M(e)) = 0. (3.10) 

Define operators p m , q m (m = 1, ■ • • , M) by 

p m = ijj++ ip m , q m = + ip m . (3.11) 

Set 

M M 
P(0) = E Vme~ m \ q(0) = 1me- ime . (3.12) 

m=l m=l 



10 



We have 



p(O) + ep(-0) = V + (e,e)-rj(e,e), (3.13) 
+ eq{-6) = r) + (e,e)+ri(9,e). (3.14) 



The following relations hold for m > 1. 

M+l M+l 

J2 p(0^ M )e imd ^ =0, £ g(0_^M)e^— = 0. (3.15) 

/i=— M ii=—M 

Therefore, p m and g m can be obtained by Fourier transformations. 

i M+l 

Pm = 2 (M + 1) Ej r l + ^M,s)-r l (e,,M,s))e ime ^, (3.16) 

i M+l 

9m = 2(M + 1) ^ (^(^M, e ) + r/(^,M,e))e ime — . (3.17) 

V > n=-M 



We define 



for p G End ^(C 2 )®^ . We call ( p ) £j n the expectation value of the operator p. 

Lemma 3.1 The expectation values of the two products ofn + (0^M,£) and n (0 ft M,£) are given 
by 

(V + (8fJ,,Mi £ )V + (8v,M))s,N <y~(# M ,M,£)f? (0v,Mi £ ))e,N 

(V (^,M,e)^ + (^,M,e)) e ,Af (v (^,M,e)»7 (^,Af,£)) e ,JV 

/ 0i(iV- lul) \ 

= 2(M + l)(^+^,_,) U . (3.19) 

^ 2 (H-7V) y 

where —M < u,v < M and 9i(x), 62 (x) are the step function 

, 1 ,x > f 1 ,x > 

l(*)H - 2 (x) = { (3.20) 

,x < ,x < 



Proposition 3.2 T/ie expectation values of the products of p m , q m are given by 

( Pi Pm )e,N ( Pi Qm )e,N \ _ I 0~l,m ~ K e,l,n 

{ qi Pm )s,n ( qi q m ) £ ,n J \ Ke,i,m 5^ m 
where I, m = 1, 2, • • • , M and 



(3.21) 



1 J Sin 2(M+1) ( l ~ Sin 2(Af+l) (* + m ) 



7T 



^ £ ,i,m = - rr— : < 1 j, r- + e- 1 77- r~ } ■ (3.22) 

M + 1 I Sm 2(M+1) (* ~ m ) 7F Sm 2(M+1) (^ + m ) 7F 
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Proof. By direct calculation, we can check the following. 

1 M+l M+l 

( P lQmh,N = (9(M ,^ 2 E E ((ri + -r l )(e,, M ,e)( V + + V )(e u , M ,e)) e , N e iW ^e i ™ e ^ 
{Z{ivi + ijj ^ = _ Mv= _ M 

i M+l M+l 

= oTmTTT E E (<W + e^,-v)z+{N - \„\) e ^,M+me^ M ) 

^ ' n=-M v=-M 



sin • 



: (l-m)ir sm^^(l + m)TT 



- -0l >m + ITF— < j 1 TTi \~~ f • ^-^j 

M + 1 [ Sm 2(M+1) (* ~ m ) 7F Sm 2(M+1) ( / + 

, 1 , X > 0, 

Here e + (x) denotes the sign function e+(x) = < We have used the relation, 

-1 , x < 0. 

^ e+(iV - H)e^ = 2 (M + 1) a Sj0 _ . ^f 17 . (3.24) 

M =-M [ S111 2(M+l) 7rS J 

□ 

We prepare some notations. Choose < mi < • • • < m n < M and < m n+ i < • • • < m<i n < M. 
Let m^ < m 2 < ■ ■ ■ < m' 2n such that = m^) (<r G -S^n)- Define the interval I^m and Im by 
I,- M = {^e Zlrn^..! + 1 < I < m 2j }, / M = h M ^h,M U • • • U/„, M . Define t m , t/ M G End ((C 2 ) 0M ) 
by i m = qiPi ■ ■ ■ qmPm, ti M = *mi ■ ■ ■ *m^ n ■ Define R £iPP i M (l,m), R £jPq i M (l,m), R £tqp i M (l,m) and 
,qqlM 

(I, m) by 



Re,ppI M (l,m) R £ , pq I M (l,m) \ 1 (j>l Pm tl M ) e ,N q m tl M ) e ,N 



(3.25) 



Re,qpI M (J, m ) Re,qqI M (^ m ) J ^ ^ M ^'^ \ (ft tl M )e,N (ft Qm tl M )e,N 

where /, m = 1, 2, ■ • • , M. Define the matrix K £j / M by (K £ j M ) jjkelM = ( qj p k } £jN . 

Lemma 3.3 The expectation value ofti M is given by 

(t lM ) £ , N = detK £jM . (3.26) 

For I, m = 1, • • • , M, the following relation holds. 

Re :P pi M (l,m) + R £ , qq i M (l,m) = 0, R £ , pq i M (l,m) + R £ , qp i M (l,m) = 0. (3.27) 

Furthermore R £tPP i M (l,m), R £tPq i M (l,m), R £tqp i M (l,m) and R em i M (l,m) have simple formulas. For 
l,m G Im, the following relations hold. 



R £ , pp i M (l,m) R £! p qlM (l,m) \ = I ft,m -(K £ j M )i,m 
R £ , qp i M (l,m) R £mlM (l,m) ) \ {K £ j M )i, m Si, m 



(3.28) 
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Proof. From Wick's theorem and ( pi p m ) £>N = fy m , we obtain ( pi p rn t/ M ) e ,jv = ( */ M )e,Nk,m- 
/,From this and Wick's theorem, we can deduce 



( tl M ) £ ,N b m ,m' = { Pm' Pm U M }s,N 

= ( Pm' Pm )s,n( tl M ) £ ,N + X! ( Qx ) £ v?v( Pm qxtl M )e,N ~ ^ { Pm' PX )s,n( Pm P\tl M )e,N 

\&Im Ag/ m 

= ^2 (Pm' q\ )e,N( Pm Qxtl M )e,N = ~( ti M ) £ ,N ^ R e,pqI M { m ^){. K e,lM)\m' 



We prepare some notations. Set 



( ■ 2(n+A0+l / /. ■ Ji/i-i- \ i — j. 

7T I sin v 2L ' TT{X — X ) Sill J — JIT — n \ X + X 



K £iN (x,x') = — . 

21/ sin 2rvr(x — x ) 

Define the integral operator K £j n,j by 



2(n+iV)+l j 



sin 2xvr(x + x') 



(K e ,N,jf)(x) = j^K £jN (x,y)f(y)dy. 



xi, 



Let us denote by det ^1 — \K £ ,n,j} an d det J 1 — \K £ ^,j 

\ x 'l i ' ' ' i x n 

minant and the n— th Fredholm minor determinant, respectively. Namely, we have 



(3.29) 



□ 



(3.30) 



(3.31) 



, the Fredholm deter- 



«(i-^)=g^ I- 1 



, %i, ■■■ ,xi 

dxi ■ ■ ■ dxiK £tN 
'J \ xi, ••• ,xi 



(3.32) 



det 1 - XK £ ^j 



Xl, 

/ 



X 



1- 



where we have used 



! X n 



g(-A) /+n 

1=0 



■■■ dx n+ i 
Jj Jj 



dx 



n+l 



x K, 



e,N 



Xi, 



x±, 



, X n , X n +\ , 



i X n -\-l 



• X n , X n _)_i, ••• , X n _|_; 



,(3.33) 



K, 



Xl, 



e,N 



Set 



Xl, • • • , X; 



(3.34) 



Re,N,j(x,x'\X) = ^2\ l ■■■ dxi • • • dxiK £i Ar(x,xi)K £i Ar(xi,x 2 ) ■ ■ • K £jN (xi,x'). (3.35) 

«=0 J J 
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Define the integral operator R £: n,.j by 

(Re,N,jf)(x) = J Re,N,j(^ (3.36) 

The resolvent kernel R e ,n,j{x, x'\X) can be characterized by the following integral equation 

(l - XK e , N ,j) (l + \Re,N,j) = 1. (3.37) 



Here we present a proof of Theorem 2.1 



Proof of Theorem 2A_ First, for simplicity, we show the n = 1 case. For s± < S2, (si,S2 € 
{1,2,- •• ,M}), we have 

( (f> Sl (f>% )e,N = ^((^ +4>s 1 ){4>f 2 +(f> S2 ))e,N = g<(^ ~ • • • W>£ + ^s 2 )) £ ,N 

= ^(-!) S2 " Sl ((9*iP«i+l)(9« 1 +lP«i+2) • • • (qs 2 -lPs 2 ))e,N (3.38) 

Applying Wick's theorem and ( pj p^ ) £j jv = <5j,fc an d ( qj qk )e,N = $j,k, we can write the above as 
a determinant. 

( 0»i #J" 2 )e,JV = J .det (( p i+ i % } £iA r) . (3.39) 
/.From (|2.17| ), Pi,n,l( x i\ x i\ £ ) = lim e-+o e iV ( <j> Sl 4>f 2 )e,i+N holds, where e = and esi -> sci, es 2 -> 



Set z/ = S2 — si. From Proposition 3.2, we obtain 



— 1 / 1 f j + 1 + s\ k s 

Pi,N,L(xi\x[\e) = -— r lim v det Sj+i,k + ~G e , - ) ) , (3.40) 

2{x 1 — x\) v-kx} l<j,k<v\ ' v \ v v; 1 



where we set 



n , /x x' x -xi fsm L ] n(y + y)(x 1 -x 1 ) sin 2L > ir{y-y){x l - x x )\ 

Ge{y, y ) = = < ; — i , , ,s, , s + £ : — — ( WT, x \ • (3-41) 

L [ sm^niy + y'){x' 1 - X!) sin ^x{y - y')[x\ - x x ) J 



We apply the following relation to the above equation. 

lim v det (-S j+1>k + -XH ( J —, -)) = -(-A)tf(0, 1) - (-A) 2 C dy 2 H I ° 

v-yoo l<j,k<u \ J 1 l> \ V V ) ) Jo \ i 



> m+ i if 1 f 1 , , TT { 0, 2/2, • •■ ,y m +i 



m! 7o Jo 



(_ A) m+i / ...i d y 2 ... d y m+1 H\ — (3.42) 



1, 2/2, ••• ,2/m+l 



Here H(y\,y2) is a continuous function, and we use 



# i y ;' " y r i = 1 < d f t < ( 3 - 43 ) 

y[, ■■■ y'n I 1 ^' k - m 
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We can write down 



Pi,7V,L(ziKk) 



7T 



7T 



•El 



2 \ m+1 i r'i 

7r / ml 



+ I -- ) — r I ■■ I dy 2 ---dy m+ iK £)N f V2) '^+1 . + 

' /a:i V x[, y 2 , ■■■ ,y m+ i 



(3.44) 



Now, we have proved n = 1 case. Next we shall prove the general case. ^From Proposition |3.2| and 
Lemma I3.3L we can deduce, 



lim ( ij } £ ,n+N = det 1 - XK e NI 



(3.45) 



/.From Lemma 3J3, we see J2iei M (Ke,i M )m,iRe,qpi M (^ m> ) = $m.m' ■ Comparing this relation to the 
relation ^1 — XK^nj^ {l + A^ £) jv,/ P ^ = 1, we can deduce the following. 



lim (—) R E qpI (mj,m k ) = A R E ,N,l p (?j, x k \X)\ 

M^oo \ Li J 



(3.46) 



for rrij ^ rrik, jjrrij — > Xj (j = 1, • • ■ Choose < mi < • • • < m n < M and < m n+ i < 
■■■ < ni2n < M. Let < m' 2 < ••• < m' 2n such that = rn a (j) ( a £ SWi)- Set m'- = 
rrij , : odd, 

^From the parity argument, we obtain 

rrij + 1 , a(j) : even. 

(*mi ' • • 0m^' ' ' ' §ml ■ ■ ■ tm' 2n )e,n+N = 0, (t^ ■ ■ ■ ■ ■ ■ • • • t m ^) £ ^ n+N = 0. (3.47) 



The expectation value (4>m"4>m" ' ' ' 4>m"4> » </> » ' " " ^ » )e n+7V can be written as Pfafnan. (See 

1 2 n m n+l m n+2 m 2n ' 

p. 967 of [?]). Furthermore, from ( |3.47| ), we can write the expectation value as a determinant. 

(4>m"4>m" ■ ■ ■ <t>mg4>m» 4>t," ' " ' Ki" )e,n+N 

1 2 n rn n + l m n+2 m 2n 



tl M )e,n+N 



(_l)|n(n-l) det 

l<j,k<n 



{t m i ■ ■ ■ (f) m ii ■ ■ ■ t m 'it m " ■ ■ ■ (j>y,i ■ ■ ■ t m ' ) £ ,n+ 



N 



tl M )e,n+N 





n 




| det ( 




' l<j,k<n v 



(3.48) 



/,From the equations fl3.45|) , fl3.46j ), ( |3.48|) , we can deduce 

2n 

(-A) n det (l - XK £>N ,i p ) det (Re^xj, x' k \X)) 



f M\ n 

M^OO \ L I 12 i m „+2 m 2n 



(3.49) 
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where j^m'j — > Xj, jjm' n+ j — >■ a?'-, (j = 1, • • • , n), when M —> oo. Using the Fredholm identity, 



(-A) n det (l - XK £ ,N,i p ) det (R £ , NJp ( Xj , x' k \\)) = det I 1 - XK £ ,n,i p 

\ I l<j,k<n ' 



we can deduce the following. 



Urn (ff 



l\ n / 2 * 
2 det 



n m n + l m n + 2 m 2n 



•^1) *^2> ' 



This complete the proof of the general case. 

Fredholm minor series in this correlation function is a finite sum because 



(3.50) 



T 



(3.51) 



□ 



K 



xi, ■■■ ,x t 



for m > 2(n + N). 



'1; "' 

To see this, define an m x M matrix Am{o\x\, ■ ■ ■ , x m ) by 

(A M (a\x l ,---,x m )) hk = e takx i, for j = l,---,m, fc = -\{M - 1) 



Using this matrix, we obtain the following 
Xi, ■■■ ,Xi 



(3.52) 



,...,-(M-l). (3.53) 



s,JV 



1' 



, x\ 



/ • 2(ra+iV)+l / /x\ 



: •••'»') •' ,l( '< ( • 1 2:< 



#1, ' ' ' ,X m ) ^2(n+iV)-l I £ 



7T 



/ / 



£l,-,£m=± 

Here A T represents the transposed matrix. ^,From elementary argument of linear algebra, we can 
see det (Am (a\ x\, • • • , x m ) A T M (f3\ x\, ■ ■ • , x' m )^j = 0, for m > M + 1. Now we have proved ( p. 52 ). 



4 Generalized fifth Painleve equation 



The purpose of this section is to give a proof of Theorem 2.3. Following H , we describe the correlation 



functions in terms of the generalization of the fifth Painleve equations, which are given by Jimbo, 

N 
L 



Miwa, Mori and Sato in the thermodynamic limit (N,L — > oo, j- = Po '■ fixed). Set 



= sin P0 «(x-x>) + e S in Po n(x + x') _ 
x — x' x + x' 
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Define the integral operators K £y j by 

(K e%J f) (x) = J j K £ (x,y)f(y) dy. 

Set 

oo „ „ 

R £>J (x,x'\X) = ^X l ■ J dxi ■ ■ ■dxiK £ (x,xi)K £ (xi,x 2 ) ■ ■ ■ K e (x u x'). 

Define the integral operators R £: j by 

(R £tJ f)(x) = J j R £>J (x,y\X)f(y)dy. 
The resolvent kernel R £: j(x, x'\\) is characterized by the following integral equation, 

(1 + XR £ ,j)(l + XK £ ,j) = 1. 



Let us denote by det — XK £) j^j and det ^1 — XK £y j 
and the n— th Fredholm minor determinant, respectively. Namely, we set 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



X\, ■ ■ ■ , x n 
Xi, ■ ■ ■ , x n 



the Fredholm determinant 



<tot(l _^)=g(^/ / ../ j 



i x\, ••• ,xi 

dx\ ■ ■ ■ dx\K £ 
'J J J \ Xi, ■■■ ,Xi 



(4.6) 



det 1 - XK £ j 



where we have used 



We set 



Xi, 



X 



, X t 



> X n 



x u 



x'i, 



,xi 
, x\ 



E 

1=0 



■■ dx n+l ■ ■ ■ dx n+ i 
J j J J 



Xl, ■■■ ,X n , X n -\-\, ■■■ ,X n +l _ 

x Ke | t f 1 , (4.7) 



X 



1> 



> x n , 3^n+l) 



d. e ,t (K E ( Xj ,x' k )) . 



i<j,k<l 



(4.8) 



Pn( x li " " " j x 'n\ x 'li ' ' ' i x n\ £ ) ~ nm „ Pn,N,L(x'\i " " " j x nl x l J ' ' ' > x nl £ ) 



(4.9) 



N,L^ca, — =po 

Contrary to the case in a finite box, the Fredholm minor series in correlation functions is infinite 
series and correlation function p n (x' 1 , • • • , x' n \x", • • • , x'^\e) becomes a transcendental function. In the 
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sequel, we shall study the differential equations for correlation functions in the thermodynamic limit. 
In what follows we can choose such a scale that 7rpo = 1- 

We prepare some notations. Let — oo < a\ < a 2 < • • • < a 2m < +00. We denote by / the interval 
defined by I = [ai, a 2 ] U • • • U [a 2m _i, a 2m ]. Set 



, sin(x-x') 
L(x, x ) — 



x — x' 



Define the integral operators Lj by 



(Ljf)(x)= J i L(x,y)f(y)dy, 



Set 



Si(x,x f \\) = y^ \ l ••• / dxi • • • dxiL(x, xi)L(xx 7 X2) • • • L(xi,x f ) 

^ ji Ji 



Define the integral operator Sj by 



(S I f)(x)= / S I (x,y\X)f(y)dy. 



The resolvent kernel Sj(x,x'X) is characterized by the following integral equation, 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



[1 + XSi)(l - XL Z ) = 1. 



(4.14) 



Set 



Si 



xi, 



Xi 



x \1 ' ' ' 1 x l 



l<J,fc<( 



(4.15) 



Set 



1 f (x - ai)(x - 03) • • • (x - a 2n _i) 
2m I (x - a 2 )(x - a 4 ) • ■ • (x - a 2n ) 



Sf(x,x'|A) = f>< / ••• / dj/i-.-dy, 

;=o 1 1 



* — ^h I (yi)L(y 1 ,y 2 )hi(y 2 )---L(yi_ 1 ,yi)hi(yi)L(yi,x'), (e = ±) 



(4.16) 



(4.17) 



2i(x - y x ) 

where the integration § c dy^ is along a simple closed curve Cj oriented clockwise, which encircle the 
points ai, • • • , a 2m . In ( 4.17D , x is supposed to be outside of C/. We denote <Sj (x, x'|A) those obtained 
by letting x inside of Cj in (4.17). Si(x, x'\X) is an entire function in both variables x, x'. Sj(x, x'\X) 
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is holomorphic except for a pole at x = x'. Sj(x,x'\X) has branch points at x = a±, ■ ■ ■ ,02m- The 
singularity structure of Sj(x,x'\X) is as follows. 



Sf(x,x'\X) - Sj(x,x'\X) =enXh I (x)Si(x,x'\X). 



(4.18) 



We set 



s £j (x\\) = jr\ i ( ... / d yi ---d Vl 

, n JCr JCr 



1=0 JC I J °I 

x L(x,y 1 )h I (y 1 )L(y 1 ,y 2 )h I (y 2 )---L(y l . 1 ,y l )h I (y l )e £iy ', (e = ±), (4.19) 



and 



00 

^f/(^|A)=E A ' / ••• / dyi---d yi 

1=0 1 1 



£ i e e'i(x-yi) 

x — rh I (y 1 )L{y 1 ,y 2 )h I (y 2 )---L(y l _ l ,yi)h I (yi)e £tyi , (e,e' = ±), 



(4.20) 



2i(x - yx) 

In ( 4.20 ), x is supposed to be outside of Cj. We denote by S £ £ 'j(x\X) those obtained by letting x 
inside of Cj. The singularity structure of S £ j(x, x'\X) is as follows. 



5 e /(x|A) - S £ £ i(x\X) = EKXh^x) S £j (x\X). 



(4.21) 



We define the matrices Yj(x), Yj(x) by 
Y T (x) = 



S +I (x\X) s-jixix) - = S +I {x\X) s; 7 (*|A) 
S-t{x\X) SZMX) ) ' \ S_j(a?|A) 5l 7 (x|A) 

^From the relation ( 4.19| ), we obtain the following monodromy properties. 



(4.22) 



1A 



Y T (x) = Yj(x) 



(x - ai)(x - a 3 ) • • • (x - a 2m -i) 1 \ 



(4.23) 



(x - a 2 )(x - a 4 ) • • • (x - a 2m ) 

The matrix Yj(x) is holomorphic and detY/(x) = 1. It is known that Yj(x) satisfies the linear 
differential equation ( 4.49 ). See |§. We define the matrices Yj a ' a \x) and Yj a,a -*(x) by 



Y^ a '\x) = {x-a){x-a') 



Sj(x, a\X) Sj (x, a\X) 



Y, 



(a,a') t 



x) 



Si(x,a'\X) Sj (x,a'\X) 

Sj(x,a\X) (x — a)(x — a')Sj (x, a\X) 
Si(x,a'\X) (x — a)(x — a')Sj(x, a'\X) 



(4.24) 
(4.25) 
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/,From ( 4.17| ), we obtain the following formula. 



1 \ / §A 



The matrix y^ a ' a ) (x) is holomorphic and 



detY r (a,a ' ) (s) = ^— A 5/ (a, a' | A). (4.27) 

Define the matrix (x) by 

v(o,a')/ n v(a,a>), J 1 °\ f w J S+(x,a|A) S7(x,a|A) . 

00 W = F ' W i-1 1 ] =(I ~" )(X ~ 0) U + MA) Sr(x,a'|A) I" (428) 
The matrix y/ 1 ^ ^(x) has the following local expansion at x = 00. 



where 



iffW = (sftf + O Q)) x exp , * ° , (4.29) 



«s*> = ( s -' (a|A) S+ ' ( "! A> I ( : 4 ] . (4. 3 o) 



5_/(a'|A) 5+/ (a' I A) / V 



2 



We set 



Zj a ' a,) (x) = 5#£' ) Y/^x), Z} a ' a,) (x) = 5}^ V^^x). (4.31) 
4*' a ' } (x) is so normalized that the local expansion at x = 00 takes the form 

Zjgp{x) = (l + O Q)) x exp L ^ M j . (4.32) 

Here we start to consider our problem for correlation functions. Let < x' x < ■ ■ ■ < x' n < 00, < 
x'l < ■ ■ ■ < x'n < 00. Let I p the union of n intervals I p = [xi, xq\ U • • • U [x2 n -i, X2 n ], where < xi < 
• • • < X2n < 00 is the re-ordering of x' l5 • • • , x' n ,x'(, ■ ■ ■ , x" n . Set I n = [—X2 n , — X2n-i] U • • • U [— X2, — xi]. 
In the sequel, we consider the case m = 2n, at = —X2 n , • • • , Q2n = — xi, £i2n+i = xi, • • • , a^ n = X2 n - 
We set I = IpUln- 
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Lemma 4.1 The resolvent kernel has the following symmetries. 

s! pUIn (x,-x'\\) = sj-; UIn (-x,x'\\), sf pUIn (x,-x'\\) = Sj; UIn (-x,x'\\), (4.33) 

S £ ,i p uiA-x\X) = S. £jpUln (x\X), (4.34) 

<W"*I A ) = S -e,i P uiM X )> S![ IpUln (-x\X) = SZl' JpUln (x\X). (4.35) 
The following is the key lemma. 

Lemma 4.2 The resolvent kernel has the following linear relation. 

R £jp (x,x'\X) = S IpUln (x,x'\X)+eS IpUln (x,-x'\X), (e = ±). (4.36) 

Proof. The following characteristic relation holds, 

S IpUln (x,x'\X) + X ( {S IpUln (x,y\X)L(y,x>) + S IpUln (x, -y\X)L(-y,x')}dy = L(x,x'). (4.37) 

J Ip 

^From ([4.37] ) and the relation e 1 = 1, we derive the following characteristic relation, 

S IpUln (x,x'\X) + eS IpUln (x,-x'\X) (4.38) 

+ A / (S IpUln (x,y\X) + eS IpUln (x,-y\X))(L(y,x') + sL(y, -x'))dy = L(x,x') + eL(x,-x'). 
J ip 

This means the equation ( [4.36| ), □ 
Let us derive a formula for d logdet ^1 — XK £ j p s j . 

Proposition 4.3 We set ui £ j p (X) = d logdet (1 — XK £ j\ . Then we have 



2n d 
^e,i P W = trace ( J2 YpUiAxj)' 1 ^- Yj p \ji n ( 



Xj 



dxj 



X=Xj \ 

e ^tram\J2Yi p ui n (x j )- 1 Y IpUIn (-x j )\ M (4.39) 



2n 



trace I /J <5AjAfc^4(xj)A(<5xfc) d log(xj — fefc) 

, ,5=± l<j<fc<2n 



2/i 



/fere i/ie matrix A(xj) is defined by 



e trace ( ^ XjA(xj) \ A^dxj + \ \ ^ J ^ ~ A i^(- X i)^ J" ] ■ ( 4 - 40 ) 



Afo) = Y IpUln ( Xj ) I ° 1 ) i^uijx,-)- 1 , A, = (-iy +1 ^. (4.41) 
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Proof. It is easy to see that 



log det (l- Ait A 



= {-iy +1 \R £ j v (x 31 x 3 \\) 

= (-l)- ?+1 A{5 /pU/n (a; i ,a; i ) + £S IpU i n (xj, -Xj)} . 
^From the definition, we can derive the following formula. 

S+i(x) S+i(x') 



det 



S-Ax) S-Ax') 



2i(x — x')Si(x, x). 



Using this formula, we obtain 

Si(x,x) 



8 S +I (x) 



idetf S+l(X > * 
1 1 S-i{x) &S-i(x 



% I ~ _■, d ~ I 1 
-trace Y I (x) —Yi(x) 
2 \ dx \ 



1 / S+i(x) S +I (-x) 
bi(x,—x) = - — det 

\ S-^x) S-i(-x) 

I / ~ , - / 1 

—trace YAx^YA-x) 
Aix \ \00 



Hence we have the first line ( |4.39| ). Substituting ( [4.22 ) into the differential equation, 



2m 



dYi(x)Yj(x) 1 = 2_. ^jA(aj) d log(x — aj) + A^dx, 

3=1 



(4.42) 
(4.43) 



(4.44) 



(4.45) 
(4.46) 
(4.47) 
(4.48) 



(4.49) 



which was derived in and comparing the coefficients of dx at x = Xj, we obtain the following. 







Xj £= ± fc=1 X j fcX fc 



->/ P u/ n (^) < 



2tj 



E E ^ . £ _ A ^, L fc " } ^u/J^), 



e=± fc=i X 3 EXk 



(4.50) 



where A r 



i 
-t 



Substituting this relation into the first line ( [4.39 ), we obtain the second 



line Q4.40| ). □ 
Remark. It is known that the matrices A(xj) are solutions of the generalized fifth Painleve 
equations introduced in S. 
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y 
„./ 



Let -oo <y,y' < +00 ,(y ^ y'). In 



Let us derive a formula for d logdet ^1 — \K £ j p 
the sequel, we distinguish the following four cases. 

i- y, y' + ■ ,X 2n - 

2. y' ^ xi, ■ ■ ■ , x 2n , y = Xj for some j. 

3. y^xi,---, X2n,y' = Xj> for some j' . 

4. y = Xj,y' = Xj' for some distinct j,j'. 
Set 

1. J(y, y') = {0,1, ■■■,2n + l}. 

2. J(y,y') = {0,1, ■■■,2n+l}\{j}. 

3. J{y, y') = {0,1, ■••,2n+l}\{/}. 

4. J(y,y / ) = {0,l,---,2n + l}\{i,j / }. 
Here we set xo = y, X2n+i = y'. Set K(y, y') = J(y, y') \ {0, 2n + 1}. Set the notations as follows. 

M;, ° ' • (i = ±l,---,±2n), M = M 2n+1 = ' " ■ ( 1.51) 

\ / 



A (,y) = (_i)i+iAl(y - x,)G/ - x,-), A^f' } = (-1)' A^y + x.Xy' + x,) , (j = 1, • • • , 2n), 
A? y) = = 1. (4.52) 



Let us state the Proposition. 



Proposition 4.4 We set oj { ^\\) = d logdet I 1 - XK £jp 



y 
y' 



We denote by d the exterior 



differentiation with respect to Xj (j G J(y,y')). Then we have 

= J2 S (l + A(^)) 1 (rtf y,) - e^- &y,) -dy- dy>) . (4.53) 
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Here we set 



n 



(y,y') 



trace[ £ ^ Z^)^ ^£(x) 



Mjdxj 



( 



trace 



J2 B^B^ dhgixi-xj) 



E E ^r^? 



^ 4" 



dx, + £ >?/ MoodXi ] , 

ieJ(y,y') 



(y,y') 



trace [ £ (^"^(-^M^ ) , 



_/±(y-y') 



(4.54) 



(4.55) 



(4.56) 



(4.57) 



i-y,y') n(-y,y') Q (-y,y') - 1 ^ Q (y,-y') p,(y-y') Q (y,-y') ^\ M ^ dx j 



trace I V 5™ J S™ ^ + Sf^ »' " y 



where 



A (y,y) = det 5}; 



:(«/,-«/') _1 <? (w/') 



I p Ul„,oo J I p Ul n ,oo I ' 



(4.58) 



£(!/,✓) = Af^'^^^OM^^Jx,)" 1 , (j = ±l,---±2n, 0,2n + l). (4.59) 



Proof It is easy to see that 



d^- 



log det l-Ai^ Mp 



y 
y' 



(-iy +1 x 



Re,l p 


( y xj 






V V' x i 





i? £ ,i p (y,y'|A) 

Then lemma L2 and the following imply the jth part of fl4,54 ). 

(y - ^Ky 7 - x k ) 



, (j^0,2n + l). (4.60) 



SlpUl„ 





Xj 











(y - y'X^j - ^fe) 



SlpUl„ 



y y 



For j = 0, 2ro + 1, the following imply the j — th part of (|4.56|). 



dy 



logdet 1-A^ £)7p 



y 



<9y 



logi? £iJj ,(y,y'|A). 



(4.61) 



(4.62) 



The second lines ( [4.55| ), ( |4.58| ) follows from the first ones by the same argument as in Theorem 4.3. 

□ 
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(v v') (—vv 1 ) (—vv r ) (—vv') ^ 

Remarks. It is known that the matrices B) and S) , , T L-.B- S T , ,' T ' are solutions 

of the generalized fifth Painleve equations in ||. For special cases y = Xi, y' = Xj, we have the 

following formula A^- V,y '^; 

x i ~ x i Ajxuzi) = trace (A(xi)A(- Xj )) Q ^ 
Xi + Xj 1(01 



trace | A(xi)A(xj) 

Finally we give a proof of Theorem [2,3j . 
Proof of Theorem |2.3| Use the following formula 



1 



R 



e,I„ 



yi ■■■ Vk 
y'i ••• y'k 





yi ■■■ 


yk \ 


det |l - XK eJp 




y[ ■■■ 


y'k ) 


(-A) fc det(l-AK £ , /p ) 



(4.64) 



and apply Proposition 4.3 and Proposition |4.4j . □ 
For n = 1 and = x' < x case, because -R[o jX ](0, x|A) = 25[_ XjX ](0, x\X), the differential equation 
becomes simpler form. 

d d 
— logpi(0|x|+) = — logS^^ylA) 

= trace ^50(0, -x.x) + ^i(0,-j;,j;)^ + ^001^0(0, - 1. (4.65) 
( i \ 

Here = . The 2x2 matrixes Bj = jB,-(ooj oij 02), (i = 0,1,2) depend on three 

parameters ao, 0,1,0,2 and satisfy the following differential systems that have the singularities at 
y = ao, ai, 02, 00. We denote by d the exterior differentiation with respect to y, ao, a\, a®. 

dZ [aua*](y) = Wlog(y - a ) + B 1 dlog(y - 01) + 5 2 dlog(y - a 2 ) + A^dy) Z { { ^(y), (4.66) 

where the 2x2 matrices Z^'^ (y) are defined in ( f4.31| ). The integrability condition 

d(d Z\ ai ' a2 hy) Z^iiy)- 1 ) = dz\ ai ' a2 hy) Z^iy)' 1 A z\ aua2) } (y) Z^hyY 1 , (4.67) 

V [aa,a2\ ^ y ' [ao,a'2\ yt> 1 ) [00,02] Ky ' [00,02] va ' [00,02] va 1 [00,02] vw ' ' v ' 

gives rise to the following closed differential equation. 

2 

dB i = -J2[B i ,B j )dlog{a i -a j )-[B i ,A 00 )da i , = 0,1,2). (4.68) 

3=0 
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The eigenvalues of Bq, Bi is (0, 1). The eigenvalues of B\ is (0, 0). The diagonal of Bq + B\ + B2 
is (1,1). ^From the above matrix properties, we reduce ( 4.68|) to the Hamiltonian equations ( |2.30| ), 
( |2.31| ) and ( 2.32| ) which was introduced in ||. And we have the equation ( [2.291) , 
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